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A six-quark state with the benzene-like color structure based on a color string model is proposed 
and studied. Calculation with quadratic confinement with multi-string junctions shows that such a 
state has a ground state energy similar to that of other hidden color six-quark states proposed so 
far. Its possible effect on NN scattering is discussed. 
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PACS numbers: 14.20.Pt, 12.40.- 



In the QED world, there are almost countless struc- 
tures of matter. Among organic molecules, there are very 
interesting varieties such as methane, benzene, fullerene, 
etc. In the QCD world, up to now, only very lim- 
ited structures of matter have been discovered: quark- 
antiquark mesons, three-quark baryons, nuclei and neu- 
tron stars. However, the QCD interaction is richer and 
more varied than QED so one would expect there to 
be even more varieties of QCD matter. In pentaquark 
studies, various color structures have been proposed: 
color singlet hadron molecules [K(qq)N(g r3 )], color anti- 
triplet diquarks [(qq) (qq)q\, quark methane [q 4 q], etc. Al- 
though the high statistics experimental data of JLab and 
COSY [l| did not reveal a pentaquark signal, the the- 
oretical view point of multi-quark structures has been 
broadened in the course of pentaquark studies. 

Lattice QCD calculation does not rule out the ex- 
istence of glueballs, quark-gluon hybrids, multi-quark 
states, etc. An interesting question is whether there is 
a peculiar property of low energy non-perturbative QCD 
that is responsible for hiding or even excluding these oth- 
erwise possible structures of matter. Lattice QCD calcu- 
lations for mesons, baryons, tetraquarks and pentaquarks 
reveal a flux-tube or string like structure H, ||. The 
naive flux-tube or string model 0,0], after adjusting the 
model parameters, gives qualitatively similar estimates of 
ground state energies of hadrons compared to other quark 
models. We therefore use this model to study the possible 
structures of multi-quark systems. Other models, such as 
the bagrnodel quark compound model Q, potential 
model Q, Skyrmion model @, chiral soliton model 
and others [llj have all been used in the study of multi- 
quark systems. In general they give similar results and 
mutually support the existence of multi-quark systems. 

In this letter, we propose and study a new color string 
structure, a benzene-like structure, for the six-quark sys- 
tem. We estimate the ground state energy of this system, 
which we call QCD benzene, using a non-relativistic color 
string model [5[ with harmonic confinement. 

Strings are actually dynamical variables to which ki- 
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netic and potential energies can be attributed [4]. The 
potential energy is proportional to the total length of 
strings. String tensions can be estimated from the bag 
model or deduced from lattice QCD calculations. They 
can also be determined empirically by fitting hadron 
masses. We take the last, i.e. a purely phenomenologi- 
cal approach. Since we are only interested in qualitative 
properties of QCD benzene, we assume a quadratic con- 
finement potential for the color string rather than linear 
as this simplifies the calculation significantly. 

A comparative study showed the inaccuracy of this 
choice is quite small @- There are two reasons to ex- 
pect this: One is that the spatial variations in separation 
of the quarks (lengths of the string) do not differ sig- 
nificantly, so the difference between the two functional 
forms is small and can be absorbed in the adjustable pa- 
rameter, the stiffness. The second is that we are using a 
non-relativistic description of the dynamics and, as was 
shown long ago [T2], an interaction energy that varies 
linearly with separation between fermions in a relativis- 
tic, first order differential dynamics has a wide region in 
which a harmonic approximation is valid for the second 
order (Feynman-Gell-Mann) reduction of the equations 
of motion. A numerical check (see Table 1 below) con- 
firmed this expectation. 

For the six-quark system, there are four possible string 
structures as shown in Fig.l. The first three were dis- 
cussed in many papers in the past but we have not found 
any previous discussion for the last one. In these fig- 
ures, represents the position coordinate of the quark 
qi which is denoted by a black dot, y» represents a junc- 
tion where three flux tubes meet. A thin line connecting 
a quark and a junction represents a fundamental, i.e. 
color triplet, representation and a thick line connecting 
two junctions is for a color sextet, octet or others, namely 
a compound string. The different types of string may 
have differing stiffness [13 , [l4j . An inverted line repre- 
sents a conjugate SU(3) color representation. Both the 
overall color singlet nature of a quark system and the 
SU(3) color coupling rule at each junction must be satis- 
fied. Then various types of color coupling are allowed in 
QCD, including that shown in figure Id of a benzene-like 
form. 
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r 3 r 5 
b. Two color-octet baryons state 




c. Three di-quarks state 




r 4 

d. Benzene-like state 
Fig. 1. Four possible string structures of six-quark 
states. 

In the string model with quadratic confinement, the 
potential energy of the benzene-like string configuration 



for a six-quark state has the following form, 

u = \ k \ I> - yrf + K E'(y* - ytf ] - « 

where the ^ means the summation is over the adjacent 
pairs. The string stiffness constant of an elementary or 
color triplet string is k, while kn is the compound string 
stiffness. In fig. Id, each quark is connected by a different 
compound string and each string with a different dimen- 
sion of the color SU (3) g roup representation will have a 
different stiffness [lj, [l4j. However, due to quarks being 
fermions, the Hamiltonian for a multi-quark system must 
be properly antisymmetrized under quark exchange. In 
QED benzene, each carbon is connected to the neigh- 
boring carbons with a single and a double bond, but it 
has been proven that QED benzene has D^h symmetry. 
Therefore, it is natural to use a common stiffness, kK, for 
the compound string for each color configuration, after 
taking into account the symmetry requirement. 

For given quark positions rj, i = 1, • • • , 6, in the center- 
of-mass (CM) system, we can fix the position coordinates 
of those junctions y;, i = l,---,6 by minimizing the 
energy of the system. The result is, 

_ (k + 4/t 2 )(r 2 + r 6 ) + « 2 (r 3 + r 5 ) + (1 + 3k) 2 i-i 
yi ~ (1 + k)(1 + 3«)(1 + 4k) 

(2) 

and other configurations have a similar form. Thus we 
can define a set of canonical coordinates, 



Ri = 


i(r 6 + r 3 - ri -r 4 ,, 




R 2 = 


i(r 6 -r 3+ri -r 4 ,, 




R 3 = 


y^(2r 2 - r 6 - n + 2r 5 - 


r3 - r 4 ), 


R4 = 


yj(2ra - r a + r x - 2r 5 + 


i"3 - r 4 ), 


R5 = 


yi(r 2 + r 4 + r 6 - n - r 3 


~ r 5), 


R6 = 


J^ir-L + r 2 + r 3 + r 4 + r 5 


+ r 6 ). 



(3) 



Finally, the potential energy U can be written as, 

_ k ( 3/C _ o Hi i-) 3 K. .-) 

U = - R? + R 2 . + Ro 

2V1 + 3k 1 1 + K 2 1 + 3k 3 

+-1 R4 + — ;-R| 4 

1 + k 4 1 + 4k b J w 

Using this string potential in a non-relativistic Hamilto- 
nian with effective quark mass m Q , we find the ground- 
state energy of QCD benzene to be 

M 6 = 6m + ~w [ 2,/—^— + 2\ — ^— + \ L K . | , 
2 \V 1 + K V1 + 3k V1 + 4k/ 

(5) 
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where u> is equal to sjkjm. 

Furthermore, by calculating the distance between two 
adjacent quarks and the distances between quarks and 
the CM, the spatial structure of QCD benzene can be ob- 
tained. The distances between any two adjacent quarks 
can be shown to be equal and so are the distances be- 
tween any quark and the CM. These are 

((r 4 -r,) 2 ) = / (V(l + «)(l + 4«) 

2mujy/K{l + 4k) v 

+1 + 4k + V3(1 + 3k)(1 + 4k)) , (6) 
where j = i + 1 (mod 6), and 

^ = Ta — h\ ; a \ (iV(i + «)(i + 4«) 

2Amu^n{l + 4k) V 

+3 + 12k + 4v/3(1 + 3k)(1 + 4k)) . (7) 

Clearly these two distances generally are not equal, so the 
QCD benzene cannot be planar. However, we can show 
that the two equilateral triangles, namely A135 and A 246 , 
lie on two parallel planes. The length of the side of the 
equilateral triangle is 

and the distance between the two parallel planes is 

((R 135 - R 246 ) 2 ) = ^~\f^^ (9) 
zmu) V k 

R135 and R246 represent the centers-of-mass of qiqsqs 
and g294<?6, respectively. The spatial structure of QCD 
benzene is shown in Fig. 2. Projecting the A 135 plane 
onto the A246 plane, we obtain a planar benzene. 
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Fig. 2. The spatial structure of the benzene-like state. 

To check this result, the same model has been ap- 
plied to the tetraquark and pentaquark systems. Three- 
dimensional configurations are also favored as shown in 
Figs. 3 and 4, where a solid-dot and a hollow-dot represent 
a quark and antiquark, respectively, while R12 = ri ^ r2 
and R 34 = T3 \ ri represent the midpoints of ri, r 2 and 
r 3 , r 4 , respectively. Please note Figs. 3 and 4 only express 
the coordinates of the quarks; the string configurations 
are not shown. 



r. 




Fig. 3. The spatial structure of the tetraquark state. 

In Figs. 3 and 4, r 12 = — r 2 is perpendicular to r 34 
r . — r4, and ri 2 , 3 4 = R12 — R34 is perpendicular to ri 2 
and r 3 4. The lengths of ri 2 and r 3 4 can be shown to be 
equal, 

<(ri 2 ) 2 } = <(r 34 ) 2 ) = ^ (10) 

For the tetraquark and pentaquark states, the distances 
between Ri 2 and R 3 4 are respectively: 




Fig. 4. The spatial structure of the pentaquark state. 



<(Rl2-R34) 2 > = — ^±±1 (11) 

moi V k 

((R 12 -R34) 2 > = — (12) 
raw V k 

In the pentaquark state, the distance between the anti- 
quark and the CM of the 4-quark subset is 

<(r 5 ,i234) 2 > = — (13) 
mw V ok 

Consequently, multiquark states in our quark model must 
form three-dimensional configurations which is due to the 
dynamics of the systems: The string shrinks the distance 
between any two connected quarks to as short a distance 
as possible to minimize the potential energy, while the 
kinetic motion expands the distance between any two 
quarks to as long a distance as possible to minimize the 
kinetic energy: The stereo structures meet this require- 
ment better than a planar one does. 

For pentaquark systems, lattice QCD obtains almost 
degenerate energies for both planar and 3-dimcnsional 
configurations. However, the entropy of a 5q system 
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is found to be larger in a three-dimensional configura- 
tion For tetraquark systems, a three-dimensional 
tetrahedral structure is rather stable against the tran- 
sition into two mesons (l5j . 

The non-zero value of (y^^sa) 2 is due to the fact 
that the anti-quark oscillates around the CM of the four 
quarks and the value of (r5 i i234) 2 is calculated with a 
harmonic oscillator wave function. 

For comparison, we also present results for two other 
color structures of the six quark system || , 



M: 



:>,:>, 



6m 



+ 



3k 



2 + 3k 



k{2k + 0.6277) 
2k 2 + 7k + 2 



2k(k + 3.186) 
2k 2 + 7 k + 2 



M 222 = 6m + ~u ( 3 + 2 J 



(14) 



(15) 



Here, M33 and M222 are ground-state energies of the two 
color-octet, hidden color state (Fig. lb) and the three di- 
quark state (Fig.lc) respectively. 

In our calculation, there are three free parameters: The 
non-strange quark mass m and the harmonic oscillator 
stiffnesses, k and nk. The first two can be fixed from the 
masses of the non-strange hadrons. We take M3 as the 
average mass of N + A and take M2 as the average mass 
of light, non-strange mesons 

M 3 = -(iV + A) = 3?n + 3a; , (16) 

— 1/3 1 \ 3 /3 1 \ 
M 2 = - I -n + -r> ) + - [ -p+ -uj ) 

4 \4 4 7 4 V4 4 J 

3 , , 

= 2m+ -uj. (17) 

From this, we obtain the effective quark mass 

m = 2M 2 - M 3 = 0.19GeV. (18) 

The constant k of an elementary string can be eliminated 
in favor of the experimental baryon mass M3. The inter- 
junction string parameter k depends on the color dimen- 
sion, d, of the string. In the MIT bag model [HI], for 
example, one finds effectively 



«d = 



(19) 



where Cd is the eigenvalue of the Casimir operator as- 
sociated with the SU (3) color representation d on either 
end of the string. However lattice QCD calculations find 
that the stiffness, Kd, of higher dimensional color charge is 
scaled by Cd instead of \J~Cd~ The compound strings 
of six-quark systems have five possible color structures: 
3, 3*, 6, 6* and 8, so each geometric configuration may 
have different states when the color dimension depen- 
dence of the string stiffness is included. This means that 



the minimum string energy also varies in general from one 
state to another. We shall not include all these (small) 
variations, but use instead an average compound string 
constant nk. 

Numerical results for these three hidden color struc- 
tures are shown in the Table 1. To check the approx- 
imation introduced by quadratic confinement the lin- 
ear confinement results with the same multi-quark wave 
functions obtained from quadratic confinement model are 
listed in Table 1 as a comparison. 

Table 1. The masses of the six-quark systems (GeV) 



confinement 


K 


Three-di-quark 


Color-octet 


Benzene- like 


quadratic 


1.0 


2.22 


2.21 


2.19 


1.5 


2.26 


2.26 


2.25 


linear 


1.0 


2.31 


2.31 


2.38 


1.5 


2.31 


2.31 


2.37 



From the Table, it can be seen that while the differ- 
ent color structures have similar masses (around 2.2 to 
2.4 GeV), the benzene-like structure has the lowest mass 
whereas the three di-quark structure has the highest mass 
for quadratic confinement. On the other hand for linear 
confinement, the benzene-like structure has the highest 
mass, whereas the other two have smaller masses. After 
taking into account the model uncertainties, we can only 
conclude that the three hidden color states have similar 
masses. With an increase of k, the differences between 
the masses decreases. In the limit k — > 00, they have the 
same mass due to having the same flux-tube structure in 
which the compound flux-tube shrinks to zero, leaving a 
hub and spokes configuration with one junction. 

The color magnetic interaction (CMI), which is impor- 
tant for hadron spectroscopy, has not yet been included 
in our model calculation. The main reason for this is that 
we have not found a method to calculate the matrix ele- 
ment of CMI for benzene-like structure. In our previous 
pentaquark calculation, we found that the tetrahedron 
has a lower mass than the Jaffe-Wilczek di-quark struc- 
ture even after taking into account of CMI [16[ . To give 
a estimate of the effect of CMI, we calculated CMI for 
hidden color states consisted of two color-octet baryons. 
the results are in the range of 20-120 MeV for two color- 
octet nucleons or A's systems. We expect the benzene- 
like structure should have energy around 2.21-2.50 GeV, 
at 0.33-0.62 GeV above the two-nucleon mass. 

The true six quark state will, of course, involve a mix- 
ture of these differing string structures. In general, chan- 
nel coupling will reduce the ground state energy of the 
system. Therefore, for channel coupling calculations in 
QCD models, these color channels should be included. 
Unfortunately, at present we do not have any reliable in- 
formation about the transition interaction between dif- 
ferent color structures, especially from a color singlet 
hadron state to genuine hidden color multi-quark states. 
The transition interaction is very much needed for the 
study of multi-quark states. 

In nucleon-nucleon (NN) scattering, QCD benzene 
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should be one of the possible intermediate states as well 
as the other color structures shown in Fig.l. In order to 
simplify our argument, we neglect the 6 and 6* strings 
temporarily. Then, if we cut the two color 8 strings, the 
QCD benzene will shrink into two color singlet baryons 
(which we take to be two nucleons for convenience in this 
discussion). Except for exchanges in quark labels, these 
two nucleons are the initial ones, as shown in Fig. la, 
in the NN scattering process. It should be possible for 
two nucleons to transform into QCD benzene and other 
genuine hidden color six quark states due to flux tube 
fluctuations when they are close enough to be within the 
range of confinement (~ Ifm). For example, the cre- 
ation of two gluons can simultaneously excite the two 
color singlet nucleons into two color octet baryons which 
are coupled together to be overall color singlet - this is 
nothing else but a QCD benzene. If the scattering en- 
ergy is near that of these hidden color states (2.2 to 2.5 
GeV), then the intermediate six quark state should be 
dominantly in a hidden color state because the number 
of hidden color states is much larger than color singlet 
one so the entropy of hidden color state is much higher 
than the color singlet one and it cannot decay into two 
colorful nucleons directly due to color confinement. It 
must transform back into two color singlet nucleons be- 
fore decaying. Such a process is similar to compound nu- 
cleus formation and therefore should induce a resonance 
around the energies of these hidden color states in NN 
scattering. We call this state a "color confined, multi- 



quark resonance" and it is different from all of the mi- 
croscopic resonances discussed by S. Weinberg [l7]. Our 
present understanding of QCD does not obviate such a 
back and forth transition between color singlet hadron 
and genuine hidden color multi-quark states. An open 
question remains as to whether or not experiments have 
observed such a color confined, multi-quark resonance. 
In fact there is a structure in the pp scattering cross sec- 
tion around 2.2 to 2.4 GeV which has been explained 
in terms of A production [l8j |. However in Our quark 
model calculations, we found that there are several iV— A, 
A — A broad resonances in this energy range, in addition 
to the benzene-like structure. Can it also be described 
by broad overlapping resonances due to these structures? 
This possibility seems not to have been ruled out and it 
is under study in our group. 

The possible existence of benzene-like structure is 
based on the description of quarks interaction via flux 
tubes. Although the lattice QCD calculation gave some 
supports to the idea, its validity is not beyond doubt. 
Other studies are needed to check if the QCD benzene 
structure obtained from color flux tube model is really a 
QCD result. 
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